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ABSTRACT 

We examine the viscosity associated with the shear stress exerted by ions in the presence of a tangled 
magnetic field. As an application, we consider the effect of this mechanism on the structure of black hole 
accretion disks. We do not attempt to include a self-consistent description of the magnetic field. Instead, 
we assume the existence of a tangled field with coherence length i cob , which is the average distance 
between the magnetic “kinks” that scatter the particles. For simplicity, we assume that the field is self- 
similar, and take A coh to be a fixed fraction f of the local disk height H. Ion viscosity in the presence of 
magnetic fields is generally taken to be the cross-field viscosity, wherein the effective mean free path is 
the ion Larmor radius X L , which is much less than the ion-ion Coulomb mean free path X u in hot accre- 
tion disks. However, we arrive at a formulation for a “hybrid” viscosity in which the tangled magnetic 
field acts as an intermediary in the transfer of momentum between different layers in the shear flow. The 
hybrid viscosity greatly exceeds the standard cross-field viscosity when U/AJ > where A = 

+ )”* is the effective mean free path for the ions. This inequality is well satisfied in hot accretion 

disks, which suggests that the ions may play a much larger role in the momentum transfer process in the 
presence of magnetic fields than was previously thought. The effect of the hybrid viscosity on the struc- 
ture of a steady-state, two-temperature, quasi-Keplerian accretion disk is analyzed. The hybrid viscosity 
is influenced by the degree to which the magnetic field is tangled (represented by = A C0 JH) y and also 
by the relative accretion rate M/M E , where M E = L^Jc 2 and L E is the Eddington luminosity. We find 
that ion viscosity in the presence of magnetic fields (hybrid viscosity) can dominate over conventional 
magnetic viscosity for fields that are tangled on sufficiently small scales. 

Subject headings: accretion, accretion disks — black hole physics — magnetic fields — MHD — plasmas 


1. INTRODUCTION 
1.1. Background 

Viscosity in accretion disks around compact objects has 
been the subject of investigation for nearly 20 yr (for a 
review, see Pringle 1981). It was recognized very early on 
that ordinary molecular viscosity cannot produce the level 
of angular momentum transport required to provide accre- 
tion rates commensurate with the observed levels of emis- 
sion in active galaxies, quasars, and galactic black hole 
candidates (Shakura & Sunyaev 1973). Consequently, the 
actual nature of the microphysics leading to viscosity in 
such flows has been the subject of a great deal of specula- 
tion. For plane-parallel flows with shear velocity u = u(y)z y 
the shear stress is defined as the flux of z-momentum in the 
^-direction. In lieu of a detailed physical model for the 
process, the work of Shakura & Sunyaev (1973) led to the 
embodiment of all the unknown microphysics into a single 
parameter a, defined by writing the shear stress as 



for plane-parallel flow, or, in an accretion disk, 

aP=^D kep , , (1.1b) 
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where P is the total pressure, rj is the dynamic viscosity, and 
fi tep , is the local orbital frequency inside a quasi-Keplerian 
accretion disk. Note the appearance of the negative sign, 
which is required so that r\ is positive-definite. Order-of- 
magnitude arguments advanced by Shakura & Sunyaev 
(1973) lead to the general conclusion that 0 < a < 1. This 
simulated the development of a large number of theoretical 
models in which a is treated as a free parameter; in many of 
these models a is taken to be a constant. This has been 
partially motivated by the fact that in quasi-Keplerian acc- 
retion disks around black holes, observational quantities 
like the luminosity depend only weakly upon ol. This 
enabled progress to be made without precise knowledge of 
the microphysical viscosity mechanisms. However, this does 
not eliminate the need for an understanding of these mecha- 
nisms, and without such an understanding, much of the 
high temporal resolution data being collected by space 
instrumentation cannot be fully interpreted. Several pro- 
cesses have been suggested to explain the underlying micro- 
physical viscosity mechanism. Initial developments focused 
on the turbulent viscosity first proposed by Shakura & 
Sunyaev (1973), and later investigated more rigorously by 
Goldman & Wandel (1995). Although the presence of turb- 
ulence in accretion disks is probably inevitable, it is unclear 
whether this particular viscosity mechanism will dominate 
over other processes that may be operating in the same 
disk, such as radiation viscosity (Loeb & Laor 1992), mag- 
netic viscosity (Eardley & Lightman 1975), and ion viscosity 
(Paczynski 1978; Kafatos 1988). 

The paper is organized as follows. In § 1.2 we provide a 
general introduction to ion viscosity in accretion disks. In 
§ 1.3 we give a heuristic derivation of ion viscosity in the 
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absence of magnetic fields. In § 1.4 we discuss the cross-field 
ion viscosity in the presence of magnetic fields. We also 
discuss the topology of magnetic fields embedded in accre- 
tion disks in § 1.5, and conclude that tangled fields are likely 
to exist. In § 2 we derive the hybrid viscosity due to ions in 
the presence of tangled magnetic fields for the general case 
of a plane-parallel shear flow. We apply our results to 
quasi-Keplerian two-temperature accretion disks in § 3 and 
discuss the main conclusions in § 4. 


1.2. Ion Viscosity 

Ion (plasma) viscosity in accretion flows has been pre- 
viously investigated by Paczynski (1978), Kafatos (1988), 
Filho (1995), and (Katz 1991). In this process, angular 
momentum is transferred between different layers in the 
shear flow by ions that interact with each other via 
Coulomb collisions. The mean free path for the process is 
then the Coulomb mean free path. Few detailed astro- 
physical models have been constructed using the plasma 
viscosity as the primary means for angular momentum 
transport because of the presumed sensitivity of this mecha- 
nism to the presence of magnetic fields. The effect of the 
magnetic field is particularly important when the ion gyro- 
radius is less than the Coulomb mean free path and the 
orientation of the local field is perpendicular to the local 
velocity gradient, because in this case different layers in the 
shear flow cannot communicate effectively. This point was 
first raised by Paczynski (1978), who argued that even for 
very weak fields (as low as 10 “ 7 G), this effect is enough to 
almo st completely quench the ion viscosity. This is prob- 
lematic, since it is very reasonable to expect near- 
equipartition magnetic fields to be present in an accretion 
flow, with strengths many orders of magnitude greater than 
10“ 7 G. 

Implicit in Paczynski’s argument is the assumption that 
the local magnetic field is exactly perpendicular to the local 
velocity gradient. However, near-equipartition magnetic 
fields would probably be tangled over macroscopic length 
scales, as evidenced, for example, by simulations of the non- 
linear stage of the magnetic shearing instability 
(Brandenburg et al. 1995; Hawley, Gammie, & Balbus 
1995; Matsumoto & Tajima 1995). We argue below that the 
presence of tangled magnetic fields effectively eliminates 
Paczynski’s concern, because ions are able to transfer a 
significan t fraction of their momentum by traveling along 
field lines connecting two different layers in the shearing 
plasma We address the issue of tangled fields being embed- 
ded in the accreting plasma further in § 1.5. 


1.3. Field-Free Coulomb Viscosity 

Consider a field-free plasma with Coulomb mean free 
path X tt and shear velocity distribution u = u(y)z, where we 
set u(0) = 0 without loss of generality. The shear stress is 
equal to the net flux of z-momentum in the y-direction. In 
terms of the field-free dynamic viscosity rf r{ the shear stress 
is given by 


du mr / kT. du , 
' ,,lt dy = ~ N Sfe ’ 


( 1 . 2 ) 


where N, is the ion number density, m, is the ion mass, and 
T, is the ion temperature (Mihalas & Mihalas 1984). The first 
factor on the right-hand side of equation (1.2) represents the 
unidirectional particle flux crossing the v = 0 plane, and the 


second factor is the magnitude of the average z-momenty 
carried by particles originating a mean distance from t 
plane. The factor of 2 accounts for the transport of partic 
in both directions across the plane. For pure, fully ioniz 
hydrogen, we have 


Tf 


K ~ v rmi I a ~ 1 - 8 x 10 N inA ’ 


(1 


where InA is the Coulomb logarithm, ^rtnj (IkTM)" 2 
the root mean square velocity of the Maxwellian distrit 
tion, and 


t u = 11.4 


77 


3/2 


AfjlnA 


(1 


is the mean time between Coulomb collisions. This yiei 
the standard result for the field-free dynamic viscos 
obtained by Spitzer (1962), 

7 * 5/2 

r\ u = 2.2 x lO" 15 -p-r- g an' 1 s -1 . (1 


Equation (1.5) is valid provided the gas is collision 
which in this case requires that the mean free path of t. 
protons X u be much smaller than any macroscopic leng 
scale in the problem. It turns out, however, that for g. 
accreting onto a black hole, XJR and XJH can exceed uni 
in general, where R is the local radius and H is the local di 
height. In this case, regions that are separated by distanc 
larger than the characteristic length over which the veloci 
varies [ v/(dv/dR ) ~ R] can easily exchange particles ar 
therefore momentum as well. In such “ nonlocal ” situatior 
the shear stress is no longer simply proportional to the loc 
velocity gradient, and one must solve the full Boltzmar 
equation in order to study the dynamics of the flo 
Another problem that arises involves the shape of the ic 
velocity distribution. When the ions are not effectively co 
fined to a small region of the flow, the local velocity di 
tribution can become distinctly non- Maxwellian due to t! 
influence of processes occurring far away in the disk, 
such circumstances, the very existence of the ion ter 
perature must be called into question. 

If this were the whole story, then the construction of dis 
models using ion viscosity would present formidable cha 
lenges. However, so far we have completely neglected tf 
effects of the near-equipartition, tangled magnetic fie 
likely to be present in an actual accretion disk. As we argi 
below, the presence of such a field will completely alter tf 
conclusions reached above if the coherence length of tl 
field is much less than the local disk height H, because the 
the ions will be effectively confined to a region of plasrr 
with characteristic size L <€ H. If the disk is thin (H < R 
this also implies L < R. 


1.4. Cross-Field Coulomb Viscosity 

Next we consider the shear stress exerted by ions inside 
plasma containing a magnetic field oriented in the 
direction and moving with velocity u = u(y)z , whe: 
u(0) = 0. Hence the magnetic field is exactly perpendicuk 
to the local velocity gradient. In hot accretion disks, or 
generally finds that X L X u for near-equipartition magnet 
fields (Paczynski 1978), where 

jt L « 0.957y' 2 fl’ 1 (1- 
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is the Larmor radius of the ions in the presence of a mag- 
netic field B. The shear stress is therefore given by 


du 

-n l T = ~2N t 

dy 


kT du , A, 
— L m, — A l rr . 
27cm, dy L A it 


(1.7) 


where ri ± is the cross-field viscosity. This is similar to 
equation (1.2), except that the magnitude of the average 
z-momentum carried by particles crossing the plane is now 
~(du/dy)A t m t because the particles originate at a mean dis- 
tance ~ a l from the plane. Another modification is the addi- 
tion of the factor ( AJA U ) which accounts approximately for 
the efficiency of the momentum transfer process. To under- 
stand the efficiency factor, imagine an ion originating on the 
right side of the plane, and spiraling about a magnetic field 
line. During one gyration, the particle crosses from the right 
side of the plane to the left side. Since A L <t A a by assump- 
tion, the probability that the particle will experience a 
Coulomb collision with another ion before returning to the 
right side is ~ AJA i{ , Hence this factor gives the mean effi- 
ciency of the momentum transfer process. The cross-field 
viscosity can also be written as 

°- 8) 

This expression agrees with the result for this case given by 
Kaufman (1960), to within a factor of the order of unity. We 
attribute the discrepancy to the approximate nature of our 
efficiency factor ( AJA U ), which does not take several details 
like the pitch angle of the spiraling ions into account, and to 
the fact that we take the ions to be originating exactly at a 
distance A L away. 

Since (AJA^ 2 1 even for field strengths as low as 10“ 7 G, 
Paczynski (1978) concluded that the ion viscosity plays a 
negligible role in determining the disk structure unless the 
magnetic field essentially vanishes. However, Paczynski’s 
conclusion relies upon the assumption that the magnetic 
field is exactly perpendicular to the local velocity gradient. 
We do not believe that this assumption is justified when the 
magnetic field is created dynamically within the disk, rather 
than imposed from the outside. When the field is created 
dynamically in turbulent plasma, it is reasonable to assume 
that the small-scale field varies randomly in time and space. 
In such situations, the field is tangled, and it is more useful 
to consider a new, “hybrid” viscosity, where the effective 
mean free path is limited by the coherence length of the 
magnetic field. We present a derivation of the hybrid vis- 
cosity in § 2, culminating with the expression for ?j hyb in 
equation (2.14). 


1.5. Tangled Magnetic Fields 

The picture of spatially intermittent magnetic “cells” 
arising out of a balance between amplification by the 
Keplerian shear and dissipation due to reconnection was 
invoked by Eardley & Lightman (1975) in one of the earliest 
discussions of magnetic viscosity. There have been several 
revisions of that picture since then, treating the effects of 
buoyancy and turbulence on the flux tubes. A review of this 
aspect including a comprehensive reference list can be 
found in Schramkowski & Torkelsson (1996). Furthermore, 
recent simulations of the Balbus-Hawley magnetic shearing 
instability (Brandenburg et al. 1995; Hawley et al. 1995; 
Matsumoto & Tajima 1995) indicate that the magnetic field 


will be tangled over macroscopic length scales. This has 
motivated us to assume the existence of a tangled magnetic 
field, although we do not self-consistently model the gener- 
ation of the field. The models we present in this paper 
assume the field to be tangled in a self-similar manner with 
respect to disk height In particular, we take the parameter 
f = KoJH to be a constant where A eoh is the coherence 
length of the magnetic field. While this might not be satis- 
fied in a realistic scenario, the expressions we derive for the 
viscosity do accommodate the possibility of a £ which varies 
with radius and are therefore sufficiently general. It might 
be noted that treatments of the magnetic shearing insta- 
bility suggest that the most rapidly growing eigenmode of 
the instability has a spatial size that is a fixed fraction of the 
local disk height (e.g., Matsumoto & Tajima 1995). 

2. ION VISCOSITY IN THE PRESENCE OF A TANGLED 
MAGNETIC FIELD 

In § 1.4, we considered the case of a shearing plasma 
containing a magnetic field oriented in the z-direction, 
exactly perpendicular to the local velocity gradient. In an 
actual accretion disk, we do not expect this to be the case 
very often. Instead, the direction of the field is likely to be a 
random function of position on scales exceeding the corre- 
lation length of the tangled magnetic field, which arises 
from MHD turbulence. It is therefore interesting to con- 
sider the shear stress exerted by ions in the general case of a 
randomly directed field. If A h < A ti , then we expect that ions 
moving between different layers in the fluid will spiral 
tightly around the field lines, in which case two of the com- 
ponents of the ion momentum are obviously not conserved. 
On the other hand, the component of the ion momentum 
parallel to the magnetic field is conserved until the particle 
either experiences a Coulomb collision with another ion or 
encounters an irregularity in the magnetic field. Hence the 
transfer of momentum from one layer to another occurs via 
the component of the particle momentum parallel to the 
magnetic field, and in this sense the particles act like beads 
sliding along a string, in what is commonly referred to as 
the ideal MHD approximation. 

The irregularities that scatter the ions may appear as 
either stationary “kinks” or fast, short-wavelength electro- 
magnetic waves depending on the details of the turbulence. 
If the particles interact with the field primarily via wave- 
particle scattering, then the waves must be explicitly 
included as a dynamical entity in the momentum transfer 
process. In fact, the shear stress due to the waves themselves 
may dominate the situation if the wave energy density sur- 
passes that of the particles. However, such large wave 
energy densities cannot be created if the field is generated 
dynamically within the plasma, as we assume here. Further- 
more, the relatively fast-moving ions that carry momentum 
in our picture will not often encounter short-wavelength 
electromagnetic waves with sufficient amplitude to scatter 
them very strongly. Conversely, the ions will be strongly 
scattered by encounters with long-wavelength, slow-moving 
lrinlrs in the magnetic field. We therefore ignore the dynami- 
cal consequences of the fast waves, and treat the irregu- 
larities as stationary kinks. 

If the field is frozen into the plasma, then the ion momen- 
tum will ultimately be transferred to the local gas via either 
Coulomb collisions or encounters with magnetic irregu- 
larities. The probability per unit length for either type of 
interaction to occur is proportional to the reciprocal of the 



No. 2, 1996 


BLACK HOLE ACCRETION DISKS 


associated mean free path. It follows that if the two types of 
interactions are statistically uncorrelated, then the effective 
mean free path X is given by 


1 _ J_ _1_ 

1 2f| ^coh 


( 2 . 1 ) 


where 2^ is the mean distance between kinks in the field, 
which is equivalent to the coherence or correlation length. 

We will continue to focus on the case of a plane-parallel 
shear flow characterized by the velocity distribution 


u = u(y)z , 


( 2 . 2 ) 


where u(0) = 0. To eliminate unnecessary complexity, we 
will also assume that the ions are isothermal with tem- 
perature T,. This is reasonable so long as the temperature 
does not vary on scales shorter than the ion effective mean 
free path X. It will be convenient to introduce a local polar 
coordinate system (r, 9, <f>) using the standard transform- 
ation 


x = r sin 9 cos <f> , y = r sin 9 sin <j> , z = r cos 9 , 

(2.3) 

in which case the velocity v r along the r-direction is related 
to v x , v y , and v z by 

v, = v x sin 6 cos <f> + v y sin 9 sin <p + v z cos 9 . (2.4) 

Let us first consider a case with no magnetic field. Then, 
viewed from a frame comoving with the local fluid, the local 
ions have a Maxwellian velocity distribution with tem- 
perature T,. However, viewed from the rest frame of the 
fluid located at y = 0, the distributions of v x , v y , and v x for 
particles located at an arbitrary value of y are given by 



due to the presence of the shear flow, where/tyi)^ gives the 
fraction of particles with ith component of velocity between 
v t and v t + dv h and J®* J{yddv f = 1. Since v x , v y , and v z are 
independent random variables, it follows from equation 
(2.4) that the distribution of v r is given by 

*’’> - (£?)”•" {" 2Mj [ ”' ' “W 005 *>'} ' 

( 2 . 6 ) 

Next we consider the effect of “ turning on ” a magnetic field 
oriented in the r-direction specified by the angles (0, <f>). If 
the field is so strong that X L < X tl , then the ions spiral tightly 
around the field lines. However, the component of the 
velocity parallel to the field (t? r ) is completely unaffected, and 
therefore the distribution of v r is still given by equation (2.6) 
even in the presence of a magnetic field. 

We wish to compute the y-directed flux of z-momentum 
due to particles crossing the y = 0 plane from both sides 
along the field line. It may be noted that since we assume 
u(0) = 0, layers on either side of this plane will have 


oppositely directed flow velocities. Since we expect t 
2 L X u in most cases of interest, we shall adopt the “ be 
on-string” model for the particle transport and work in 
limi t XJXft -*• 0, in which case v y and v z are given by 

v y = v r sin 9 sin <f> , v z = v r cos 9 . (.' 

Hence we ignore the components of momentum perp 
dicular to the field and consider only the transport 
momentum along the field lines. For the purpose of cal 
lating the momentum flux, it is sufficient to consider p 
tides starting out at a distance 2 from the origin. It folic 
that at the starting point 

y = X sin 9 sin <f> , (- 

and therefore 


u(y) = u'(0)X sin 9 sin 4> (- 

to first order in X, where the prime denotes differentiae 
with respect to y. The y-directed flux of z-momentum due 
particles approaching the origin from both sides of 
y = 0 plane is given by 

W<t>) = 2 f ° ( mi v z )ZNi v y f(v r )dv r ] , (2. : 

J— ao 

where the first term inside the integral is the z-momentu 
carried by the particles and the second term is the y-direct 
particle flux. Then to first order in X we obtain 


P(9, <f>) = 2m, N t cos 9 sin 9 sin <j> 

rkT, (2kT\" 2 „ m , a ■ a A 

x — - — ( u (0)2 cos 9 sin 9 sin (/> 

[2m, \itmij J 


(2.1 


which gives the shear stress as a function of 9 and <f>. T 
first term on the right-hand side describes the “thenr. 
stress” due to the stochastic drifting of particles along t 
field lines, which occurs even in the absence of a veloci 
gradient The second term gives the modification due to t 
presence of the velocity gradient. Equation (2.11) vanish 
when the field is exactly perpendicular to the velocity gr 
dient (sin 9 sin <f> = 0), which agrees with equation (1.8) f« 
the cross-field viscosity in the limit XJX,, -» 0. 

Equation (2.11) for the direction-dependent stress can 1 
used to construct two-dimensional models that treat bo 
the radial and azimuthal structure of the disk. In the 
models, the direction of the local magnetic field is a rando 
function of the radial and azimuthal position on seal 
exceeding the coherence (correlation) length 2^,,. In ord 
to construct one-dimensional models, we need to averai 
equation (2.11) over all directions to obtain the mean stres; 


<P> 




#*n. 


(2.1. 


where d£i = sin 9d9d<f>, and O<0 <ji, 0 < <p < 2n. Subst 
tuting equation (2.11) into equation (2.12) and integratir 
over 9 and <f> yields for the mean (direction-averaged) hybn 
viscosity 


= _<£> 
,,brb ~ u'( 0) 



( 2 . 1 - 


Note that the “thermal stress” appearing in equation (2.1 
is symmetric and therefore it vanishes upon integration. 
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We can also write the hybrid viscosity given by equation 
(2.13) as 

2 X 

flbyb = -p: 7f f . (2.14) 

where is the standard, field-free Coulomb viscosity given 
by equation (1.5). We see that no factor describing the effi- 
ciency of the momentum transfer process appears in the 
expression for rj hyb , in contrast to the cross-field viscosity 
given by equation (1.8). This is because in the hybrid case 
particles originating on the right side of the plane and cross- 
ing over definitely deposit their momentum on the left side, 
since the two planes in question are linked by a field line 
along which the particle “ slides.” Since rj hy Jrt {{ ~ (X/X u ) and 
~ (V>U 2 > it is clear that the hybrid viscosity will 
greatly exceed the cross-field viscosity if (X/X L ) $> (XJX^ 
which is likely to be well satisfied in hot accretion disks, as 
will be seen in § 3. This suggests that the ions play a much 
larger role in the momentum transfer process in the pres- 
ence of magnetic fields than originally concluded by 
Paczynski (1978). In § 3 we use our results to analyze the 
structure of a two-temperature quasi-Keplerian accretion 
disk with unsaturated inverse-Compton cooling. 

3. APPLICATION TO TWO-TEMPERATURE ACCRETION 
DISKS 


using a specific model for the viscosity. In our case the 
variation of a is obtained by substituting our expression for 
;j hyb into equation (1.1). 

In order to close the system of equations and obtain 
solutions for the disk structure, we must also adopt a model 
for the variation of the magnetic coherence length X coh 
which appears in the definition of the effective mean free 
path X (eq. [2.1]). We assume here that the field topology 
varies in a self-similar manner with the local height if, so 
that 

JU ■ & > ( 3 * 3 ) 

where £ is a free parameter which we set equal to a constant 
for a given model 

3.1. A Two-Temperature Accretion Disk Model 
In a cylindrically symmetric accretion disk, the relevant 
component of the stress arising from the hybrid viscosity is 
given by 

®hyb P=~ tJby b R 

which is equivalent to equation (1.1). We use equation (3.4) 
to derive oe hyb from */ hyb . Equations (A7)-(A9) in Appendix A 
and equation (3.4) jointly yield the following self-consistent 
solutions for the model: 


We consider the two-temperature, steady-state model 
first proposed by Shapiro, Lightman, & Eardley (1976) and 
adopted by Eilek & Kafatos (1983). The model assumes that 
the disk is quasi-Keplerian; i.e., the azimuthal velocity v + is 
equal to the Keplerian value and the radial velocity v r <! v 
In this model the ions and electrons are coupled only via 
Coulomb collisions and the electrons with temperature T e 
are assumed to radiate their energy away via unsaturated 
inverse-Compton cooling. In this case the two-temperature 
condition 7J > T e is satisfied if 

r> (3.1) 

while in order for the ions and electrons to achieve their 
respective Maxwellian distributions, we must have 

hccr > hi > *ee > (3*2) 

where t €i , t 9€ , and t it are the time scales for electron-ion, 
electron-electron, and ion-ion Coulomb equilibration, 
respectively, and f accr is the time scale for accretion onto the 
black hole. We will use the viscosity prescription given by 
equation (2.13), and we will assume that the coupling 
between ions and electrons occurs exclusively via Coulomb 
interactions. Hence we neglect the possibility that collective 
plasma processes might result in an additional coupling 
between the ions and electrons, over and above the usual 
Coulomb coupling, which could in principle lead to a vio- 
lation of the two-temperature condition. However, 
Begelman & Chiueh (1988) considered this possibility, and 
concluded that such collective processes are not likely to 
strongly affect the thermal structure of the disk. Equations 
(A1HA6) in Appendix A list the basic structure equations 
for the two-temperature quasi-Keplerian disk model. Equa- 
tions (A7HA9) in Appendix A constitute a list of the analyti- 
cal solutions to these structure equations, which are derived 
under the assumption that 7] > T e . These solutions have an 
arbitrary a parameter built into them, which in general can 
be treated as a constant or allowed to vary with radius 


c hyb = 147.31 S^fr 1/6 /i /3 (£0 2/3 r« 1 K ' . (3.5) 

/m V ' 3 

] = 3.38 X 10 11 5~ 1/3 /, 1/6 /i /3 ( ) R ; 1/2 , (3.6) 


T. = 


1.40 x 10 9 y 


' T«,(l + O ’ 

N t = 5.70 x I0 ll 5 ll6 f! n2 f{ u6 (jf 


(3.7) 


\ 5/6 


(3.8) 


where 


and 


H (M Y /6 

— = o.i755 -1/6 /f 5/12 fy 6 \rf') R i /4 ’ ( 3 - 9 > 

'■H-ar- (3,o) 


.V tit _ M R 

M * = lM Q yr- 1 ’ M8_ 1O 8 M 0 ’ *“GM/r' 

(3.11) 

The following two equations jointly define an implicit 
algebraic equation for determining <5 as a function of R, for 
given (f, y, MJM e ) 

5 113 (M V /3 

915.50 K' (3.12) 
t.™(l +tJ-57.88195>'’/r"'»/r'»/r J ^y' # R.- ! “ . 


(3.13) 
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We will restrict our attention to 1 > { > 0, since { ^ 1 
implies that the field is strongly ordered over macroscopic 
length scales, which violates our assumption that the field is 
tangled. Once a root for <5 is determined for a given £, it is 
used in equation (3.12), and the result obtained for is 
then used in equations (3.5) — (3.9) to determine the disk 
structure. For instance, <5 ranges between 10 -6 and 10 -7 for 
the model shown in this paper. The small value of 8 is 
indicative of the fact that the mean free path is determined 
primarily by 2 coh and not by k u . In principle, one could 
compute a disk model for a given y and any combination of 
(, M„, and M s . We consider 0.001 < M/M E < 1» where 
lk E = Lgjc 1 and L E = 4k GMm p c/<j t is the Eddington lumi- 
nosity and <jj is the Thomson cross section. Note that 
M e = 0.22 MJM a . Accretion rates that are close to the 
Eddington value are more likely to be significant from the 
point of view of observations. 


3.2. Model Self-Consistency Constraints 

For the models to be self-consistent, they have to fulfill 
the following conditions. 

(i) H/R < 1. This ensures that the disk remains geometri- 
cally thin and is assumed in deriving the analytical solutions 
listed in Appendix A. 

(ii) k/H < 1. This condition ensures that the magnetic 
field lines are confined within the disk. Since £ = k/H is a 
free parameter in our model and we restrict our attention to 
f < 1, this condition is automatically satisfied. 

(iii) k/R < 1. This assures us of the validity of applying 
the fluid approximation to the plasma. Imposing 
£ = k/H < 1 and ensuring H/R < 1 results in the satisfac- 
tion of this criterion. 

(iv) TJT t > 1. Satisfaction of the two-temperature condi- 
tion is assumed in deriving the analytical solutions in 
Appendix A. 

3.3. Comparison with Other Kinds of Viscosity 

Of the different kinds of viscosity that can possibly exist 
in the accretion disk, we assume the hybrid viscosity we 
have derived here to be the dominant form. We would like 
to compute other possible forms of viscosity a posteriori, 
and compare them with the hybrid viscosity. The hydrody- 
namic turbulent viscosity used by Shakura & Sunyaev 
(1973) is based on dimensional arguments, and, according 
to Schramkowski & Torkelsson (1995), is probably less sig- 
nificant than viscosity arising from MHD turbulence, in 
which the magnetic field plays a significant role. For rela- 
tively high accretion rates, one would expect rather high 
luminosities. Consequently, the contribution of radiation 
viscosity, which is characterized by an associated a rad , 
would be appreciable. Magnetic viscosity (characterized by 
a mmg ), which arises from the stresses associated with the 
tangled magnetic field, is another important form of vis- 
cosity with which we compare the hybrid viscosity obtained 
from our calculations. Appendix B describes how a„ d is 
calculated. We describe how a mig is calculated below. In 
particular, it may be noted that we need to adopt a different 
definition of the coherence length A coh , given by equation 
(3.15), in order to quantify 

If we consider the magnetic stress to be equal to the 
magnetic pressure P B = B 2 /(87r), the a. parameter arising out 


of pure magnetic viscosity is defined by 

__B 2 
*“ ,P “ 8k ’ 


(3.1- 


where P is the total pressure. It is customary to use tfc 
following definition for in discussions of magnetic vi* 
cosity, 


where is the Alfven speed. This formula arises from 
balance between amplification of the field by shear and di< 
sipation due to reconnection. Vertical hydrostatic equ: 
librium can be expressed as 

P = frL P iPH 2 - (3-l« 

Combining equations (3.14), (3.15), and (3.16) gives 

= (ir ) ■ (11, 

In our models, £ = X co JH is a free parameter, and equatior 
(3.17) reduces to a mag = £ 2 . 


3.4. Results 

Figure 1 shows a parameter space plot for a canonicalh 
rotating Kerr black hole (a/M = 0.998). Each point in tb 
parameter space spanned by f and M/M E represents ; 
potential model. Each of these models assumes a constan 



F ig. 1. — The ({, AST/Me) parameter space for the canonical Kerr metric 
with a/M = 0.998. Each point in the parameter space represents a poten- 
tial model. The disk is assumed to occupy the range R^ < < 50, where 

is the radius of marginal stability in gravitational radii. For this range 
of radii, models with parameters lying below the solid line fulfill a hTb > 
those with parameters lying above the dotted line fulfill H/R < 1: 
those with parameters lying above the dashed line fulfill a hyb > a^. The 
condition 7] > T t is fulfilled throughout the parameter space shown. The 
shaded region is the one in which our models are self-consistent insofar as 
conditions (i)-(iv) described in § 3.2 are concerned. 
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Fig. 2. — Results for a kyb , a,,.,, and a rld for model a (see Fig. 1), with 
M/M e = 0.0025 and { - 0.03. 


value of $ = X„JH throughout the extent of the disk, which 
we take to be R m < R m < 50, where is the radius of 
marginal stability in gravitational radii. It turns out that the 
condition TJT e > 1 is satisfied for models for the entire 
parameter space shown in Figure 1. Models with param- 
eters lying to the left of the dotted line are guaranteed to 
satisfy H/R < 1 for all radii. It follows that such models 
satisfy self-consistency constraints (i)-(iv) of § 3.2. Models 
with parameters lying below the solid line satisfy a hyb > a ma( 
for all radii, while those with parameters lying above the 
dashed line satisfy a hyb > a„ d . The shaded region, therefore. 



Fig. 3. — Results for H/R, X/R, and for model a (see Fig. 1), with 
M/M e = 0.0025 and { = 0.03. 



Fig. 4. — Results for ot^, g_. t . and 2 rld for model b (see Fig. 1), with 
M/M e - 0.0375 and { =■ 0.105. 


is the one in which our models will be fully self-consistent. It 
may be emphasized that our calculations take only a hyb into 
account; at—, and a„ d are calculated a posteriori. Figure 1 
therefore suggests the following: 

(a) For sufficiently small = X^JH (approximately 
^0.1) hybrid viscosity dominates over magnetic viscosity. 

(b) For sufficiently large accretion rates, the disk 
becomes puffy ( H/R > I) and the luminosity becomes high 
enough so that radiation viscosity dominates over hybrid 
viscosity. However, our models turn out to be optically thin 
to electron scattering in general. This means that the mean 



Fig. 5. — Results for H/R, X/R, and r„ for model b (see Fig. 1), with 
Aif/Af E = 0.0375 and { = 0.105. 
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free path for photons in the vertical direction is greater than 
the height H, and the radiation viscosity we calculate is 
therefore an overestimate of the actual value. It may also be 
noted that the breakdown of the thin disk condition implies 
a breakdown of the entire quasi-Keplerian structure and the 
radial velocity can no longer be ignored in comparison 
with the azimuthal velocity. 

These statements, however, can only be taken to be indic- 
ative of general trends. It would require a self-consistent 
inclusion of all the different forms of viscosity in order to 
obtain more concrete conclusions in this regard. We are 
concerned only with demonstrating the applicability of the 
hybrid viscosity mechanism to accretion disks in this paper. 

We next ex amin e two specific models from the parameter 
space shown in Figure 1, models a and b. They serve as 
illustrations of a two-temperature accretion disk model that 
employs hybrid viscosity. Figures 2 and 3 show some physi- 
cal quantities associated with model a, with M/M E = 
0.0025 and £ = 0.03. Figure 2 shows that a byb is around 
0.0045 and is also well in excess of a rtd and i M| . Figure 3 
shows that the disk is indeed thin and that X/R 4 1. Figures 
4 and 5 show quantities associated with model b, with 
M/M e = 0.0375 and £ = 0.105. Figure 4 shows that a byb is 
around 0.017, and it is well in excess of a mi| . It exceeds a rid 
in all but the innermost regions. Figures 3 and 5 also show 
that both the models are optically thin to electron scat- 
tering. We have verified that the radiation and magnetic 
pressures are much smaller than the gas pressure for both 
these models, thus justifying the assumption of gas pressure 
dominance. The ion temperatures for both these models are 
around 10 12 K, while the electron temperatures are around 
10 9 K. 

4. DISCUSSION 

We have derived a hybrid viscosity arising from momen- 
tum deposition by ions in the presence of a tangled mag- 
netic field. This viscosity is neither the usual Coulomb 
viscosity which arises from Coulomb collisions between 
ions, nor is it pure magnetic viscosity, which is due to mag- 
netic stresses. The tangled magnetic field plays a role in 
confining the ions, which makes the viscosity mechanism a 
local process. The field also acts as an intermediary in the 
momentum transfer between ions, in situations where the 
coherence length of the field 2 cob X lh where X u is the usual 
Coulomb ion-ion mean free path. Upon application of this 


form of viscosity to a specific disk model, we observe th; 
the self-consistency requirements limit valid models to sv.l 
Eddington accretion rates. We also observe that the hybn 
viscosity dominates over magnetic viscosity for small valui 
of £. The only restriction on the magnitude of the magnet 
field in our calculations is that it be at least so large as t 
warrant the assumption of nearly zero gyroradii for tl 
ions. 

We have entirely neglected any momentum transft 
arising out of short-wavelength plasma waves in the accn 
tion flow. Since the tangled magnetic field is taken to t 
arising from plasma turbulence, the presence of such wave 
is quite plausible, and it is one aspect of the problem w 
have neglected in our calculations. In fact, ion streamin 
along magnetic field lines, which is considered to be tt 
primary means of momentum transfer in this paper, migi 
very well be the source of such waves if their speeds ai 
super-Alfvenic. One could model an ensemble of such turt 
ulent plasma waves as a collection of plasmons, assign 
number density and mass to these entities, and investigat 
their role as intermediaries in momentum transfer. A sel. 
consistent calculation of the tangled magnetic fields arisin 
as a consequence of the presence of plasma turbulenc 
could also reveal magnetic flutter; temporal variations i: 
the local magnetic field (as distinct from the large-seal 
evolution of the fields due to dynamo action) that we hav 
also neglected. 

In the present work we have adopted a time-independen 
treatment. There have been a number of investigations o 
possible disk instabilities (Shakura & Sunyaev 1976; Pirai 
1978, for instance) which consider the presence of therma 
and viscous instabilities that could break up the disk am 
cause variations in the disk luminosity. The temperature 
dependent nature of any viscosity in which ions play a par 
(like the hybrid viscosity discussed in this paper) woulc 
result in a coupling of viscous and thermal instabilities. W< 
are currently in the process of undertaking an investigatioi 
of these aspects. It may also be noted that advection effect 
can assume importance and alter the disk structure near thi 
inner edge. 

We would like to thank Dr. Ethan Vishniac and the 
anonymous referee for several useful criticisms and sugges 
tions. P.S. also acknowledges useful discussions with Dr 
UlfTorkelsson. 


APPENDIX A 


CONSTITUTIVE EQUATIONS FOR TWO-TEMPERATURE, COMPTONIZED MODEL 


The basic disk structure equations are the same as those used in the disk structure calculations of Eilek & Kafatos (1983). 
which neglect radiation pressure : 


GMm t N,H 2 fi 
P ~ R 3 

„ (GMR) ,/2 M/ 2 

aP 4a Wh ' 


(Al. 

(A2^ 


3 GMM 


h 


= 3.75 x 10 21 m i lnAN, 2 k 


(7] - T e ) 



(A3) 


8a R 3 H 
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P = N t k(T, + ty, (A4) 

(A5) 

t = (A6) 


where <r T is the Thomson scattering cross section. The Coulomb logarithm In A is taken to be 15 in our numerical calculations 
and the function = 1 + t„. It may be noted that this is different from the form for gii^J used by Eilck & Kafatos (1983). 
The factors f lt f 2 and f 3 are the relativistic correction factors appropriate to the metric under consideration. These factors for a 
Kerr black hole with a/M = 0.998 are used in Eilek & Kafatos (1983). Eilek (1980) gives plots of fufi, and f 3 for a Kerr black 
hole. The relativistic correction factors appropriate to a Schwarzschild metric can be obtained by setting a/M = 0 in the 
expressions for /,, and f 3 . In keeping with the convention used in Eilek & Kafatos (1983), we make the definitions 
M 8 s M/10 8 M©, M„ s M/1 M© yr _1 , and R m = R/(GM/c 2 ). If we assume 7J > T t , equations (A1HA6) yield the following 
analytical solutions : 

7J = 4.99 x 10 13 tt* / 2 t' 1 «' 1 Rr 3,J , (AT) 

Af 8 

T' = 1.40 X lOVr.'CfiKO]" 1 » (A8) 

N, = 4.70 x 10 1 1 f{< 2 f 2 ~ 1/2 ill 2 a 1/2 R; 3/4 • (A9) 

It may be emphasized that a is a free parameter in the above solutions. 


APPENDIX B 


DEFINITION OF RADIATION VISCOSITY 


We use a^,, the a parameter obtained from radiation viscosity, as a diagnostic in this paper. We now proceed to define the 
manner in which we compute a rad . We follow Shapiro et al. (1976) in defining the radiation energy density using 

U ttd = (F/c)g( tj. (Bl) 

If the y parameter is taken to be equal to unity, eliminating g( r ei ) between equations (A5) and (Bl), using equation (A6) yields 

(B2) 

where F is the dissipated energy density. Equation (A3) is another way of defining F/H; in fact, F has to be equal to 
(3/8 it)GMM/R 3 for the disk to be quasi-Keplerian. Equating the right-hand side of equation (A3) to that of equation (A8) and 
assuming 7J > T e yields 

U rtd = 9.565 x 10 5 N, 7; T~ SI2 . (B3) 

We next adopt the definition of radiation viscosity given by Loeb & Laor (1992), 


8 Ur . d 
27 Njffj c 


We calculate a rMd by adopting the usual definition for a, akin to equation (3.4), 


®rad F — 


-K>fr.d 


dR ' 


(B4) 


(B5) 


It may be noted that we assume the disk to be gas-pressure dominated; although we do calculate a nd as a diagnostic tool, P in 
equation (B4), which represents total pressure, does not include radiation pressure in our calculations. 

This yields 


fill = = 6-45 x l0 33 <5 -1 lnA7V 3/2 T~ 512 , 

®hyb tfhyb 


where S is defined in equation (3.10). 


(B6) 
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